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ABSTRACT 

At the onset of high-mass star formation, accreting protostars are deeply embedded in massive cores 
made of gas and dust. Their spectral energy distribution is still dominated by the cold dust and rises 
steeply from near-to far-infrared wavelengths. The young massive star-forming region IRDC 18223-3 
is a prototypical Infrared-Dark-Cloud with a compact mm continuum core that shows no protostellar 
emission below 8 /zm. However, based on outflow tracers, early star formation activity was previously 
inferred for this region. Here, we present recent Spitzer observations from the MIPSGAL survey that 
identify the central protostellar object for the first time at 24 and 70 //m. Combining the mid- to far- 
infrared data with previous mm continuum observations and the upper limits below 8 /zm, one can infer 
physical properties of the central source. At least two components with constant gas mass M and dust 
temperature T are necessary: one cold component (~15K and ^576 M Q ) that contains most of the 
mass and luminosity, and one warmer component (>51K and >0.01Mq) to explain the 24 /im data. 
The integrated luminosity of ~177L Q can be used to constrain additional parameters of the embedded 
protostar from the turbulent core accretion model for massive star formation. The data of IRDC 18223-3 
are consistent with a massive gas core harboring a low-mass protostellar seed of still less than half a 
solar mass with high accretion rates of the order 10 -4 M yr _1 . In the framework of this model, the 
embedded protostar is destined to become a massive star at the end of its formation processes. 
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1. INTRODUCTION 



The onset of massive star formation was elusive to ob- 
servational research until recently. Very young regions 
of massive star formation contain large amounts of cold 
gas and dust observable from far-infrared (FIR) to mm 
wavelengths. A possibly luminous embedded object has 
not yet formed or is obscured by the large optical depth 
of the surrounding cold dust at near-/mid-infrared wave- 
lengths (NIR/MIR) The advent of the Infrared Space 
Observatory (ISO) and the Midcourse Space Experiment 
(MSX) allowed to ident ify large numbers of such Infrared 
Dark Clouds (IRDCs) (lEgan et al.lll998t iBacmann et all 



12005 ICarey et alJl^OOOt iSimon et al.ll2006l) . However, the 
IRDCs are not a well defined class, but these clouds are 
expected to harbor various evolutionary s tages. Adopt- 
ing th e evolutionary sequence outlined in iBeuther et alj 
(2006), IRDCs should be capable to contain genuine High- 
Mass Starless Cores (HMSCs), High-Mass Cores har- 
boring accreting Low-/Intermediate-Mass Protostars, and 
the youngest High-Mass Protostellar Objects (HMPOs). 
While the first group is important to study the physical 
conditions of massive cores before the onset of star for- 
mation, the latter two stages are essential to understand 
the early evolution in massive star formation. Since the 
evolutionary timesca le to form high-mass stars is short (of 
the order 10 5 yr, e.g- lMcKee fc Tanll2002l ). and the central 
evolving protostars are deeply embedded within their natal 
cores, observational discrimination between these different 



evolutionary stages of IRDC s is a c hallenging task. 

Recently, IBeuther et al.1 ([2005) observed with the 
Plateau de Bure Interferometer (PdBI) a filamentary 
IRDC containing a compact 3.2 mm continuum peak re- 
lated to a massive core that remains undetected in the 
SPITZER GLIMPSE survey up to 8 /zm (Fig. [TJ left 
panel). Based on the GLIMPSE non-detection one could 
speculate that this may be a genuine HMSC, however, 
there are additional indicators leading in the opposite di- 
rection. The Spitzer 4.5 /zm band shows weak emission 
right at the edge of the 3.2 mm continuum core. This so- 
called green-fuzzy emission is usually attributed to shock- 
excite d H2 emission which is espe cially prominent in this 
band (|Noriega-Crespo et al.ll2004T) . Such H2 features can 
be attributed to early outflow activity. This scenario is 
supported by single-dish CO and CS spectra showing non- 
Gaussian line-wing emission indicative of molecular out- 
flows. Furthermore, NH~3(1,1) and (2,2) inversion line ob- 
servations revealed relatively high gas temperatures of the 
order 33 K which would not be expected in t he case of a 
starless core. Therefore, IBeuther et all {2005) argue that 
we are witnessing the onset of massive star formation. 
Based on the outflow signatures and the high NH3 temper- 
atures, a central protostar should have formed already, but 
it must still be in a very early evolutionary stage to remain 
undetected up to 8 /zm. The detection of at least three 
4.5 /zm outflow features indicates that even a multiple sys- 
tem may be embedded in this co re. At a temperature 
of 33 K and a distance of ^3.7 kpc (Sridhar an et al.l l2002. 



1 Max-Planck-Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany 

2 Astronomisches Recheninstitut am Zentrum fur Astronomie Heidelberg, Monchhofstr. 12-14, 69120 Heidelberg, Germany 

3 The near kine matic distance was c hosen because IRDC at >12kpc arc unlikely to be observable. The kinematic distance uncertainty is of 
the order 1 kpc l IBrand fc Blitdll99"3h . 

1 



2 



Beuther & Steinacker 



l200ffil ). iBeuther et all (|2005h calculated from the 3.2 mm 
continuum flux a mass and column density of the gas core 
of -184 M and -10 24 cm- 2 (A v - 1000). 

In the framework of the above mentioned evolutionary 
sequence, IRDC 18223-3 should be part of the High-Mass 
Cores harboring accreting Low-/Intermediate-Mass Pro- 
tostars. For such a source, it is expected that the spec- 
tral energy distribution (SED) rises sharply at mid- to far- 
infrared wavelength and should hence become detectable 
there. The newly released Spitzer legacy survey MIPS- 
GAL (MIPS Inner Galactic Plane Survey) is the ideal re- 
source to search for such MIR/FIR emission and study 
this young source in more detail. 

2. DATA 

The MIPS 24 and 70 /jm data were taken from the 
Spitzer archive of t he recently released MIPSGAL survey 
(|Carev et all 12005). Fluxes were extracted via aperture 
photometry subtracting the background emission from 
a close-by region. The Spitzer IRAC 3.5 to 8/im and 
the PdBI 3.2 mm co ntinuum data were first presented in 
IBeuther et ail (|2005l ). Furthermore, we use a 1.2 mm con- 
tinuum flu x measurement obse rved with the IRAM 30 m 
telescope ([Beuther et all [2002T) . For the 3.2 mm data, 
we use the flux within the 50% contour level (radius 
of — 10000 AU) to avoid contamination from the large- 
scale filament. Similarly, for the 1.2 mm data we use 
the peak flux measu rement within the central 11" beam. 
IBeuther et all (|2005l ) estimated for the whole large-scale 
filamentary structure that the PdBI observations suffer 
from about 25% of missing flux. However, since here we 
are not interested in the whole filament but only analyze 
the flux from the central compact peak-emission, this effect 
is considerably lower. The accuracy of the flux measure- 
ments at mm wavelength is estimated from the data to be 
correct within —15% and for the 24 and 70 /jm within 20%. 
The 3cr upper limits of the four Spitzer IRAC datasets are 
0.05 mJy at 3.6 and 4.5 /xm 0.13 mJy at 5.8 /zm and 0.15 
at 8 /xm (|Beuther et al.ll2005h . 

3. RESULTS 

Figure [T] presents overlays of the MIPS/IRAC mid-/far- 
infrared data and the 3.2 mm observations. The most 
striking result is that the 3.2 mm core, which is dark at 
least up to 8 /zm, is now detected in the 24 and 70 /mi 
bands. Previously, we could only indirectly infer from out- 
flow indicators and warm gas temperatures that the core 
likely harbors already a very young protostar, but the new 
24 and 70 /tm data clearly identify this source now for the 
first time shortward of 100 /mi. 

Combining the 3.2 and 1.2 mm data from the Rayleigh- 
Jeans part of the spectrum with the 24 and 70 /im fluxes 
on the Wien-side of the SED, this dataset allows to derive 
the physical properties of this High-Mass Core at the onset 
of massive star formation in more detail. Table [1] presents 
the fluxes from 24 /im to 3.2 mm, and Figure [5] shows the 
resulting SED. 

The spectral energy distribution was fitted with Planck 
black-body functions accounting for the wavelength- 
dependent emis sivity of the dust. Th e assumed dust com- 
position follows IDraine &: Led <| 19841) . As a first order ap- 
proach, we tried to fit the dataset with a single-component 



black-body function. While this may work for the mea- 
surements upwards of 70 /mi, the 24 /im point shows sig- 
nificant excess emission to a single-component fit, and a 
second warmer component has to be added. Figure H] 
presents a two-component fit to the data. While most 
of the flux, mass and luminosity stems from the cold gas 
and dust with an approximate temperature of 15 K, we 
find another warm component with a temperature around 
>51 K. The gas masses associated with the cold and warm 
components are —576 and >0.01M Q , respectively. Since 
on the Wien-side of the SED the emission is not optically 
thin anymore, we consider the masses and temperatures 
of the warm component as lower limits. The mass differ- 
ence of the cold componen t between the n ew fit and the 
older mass estimates from IBeuther et all ((2005) are due 
to the lower dust temperature we now derive (15 K versus 
33 K from the NH3 gas observations) and the different as- 
sumptio ns about the dust gr ain size distributions (follow- 
ing now IDraine &: Led 1984 wher eas previously we used 
the approach by lHildebrandlll983l ). Furthermore, the in- 
tegrated luminosity one derives from this two-component 
fit is — 177L Q (174 and 3Lq in the cold and warm com- 
ponent, respectively). 

4. DISCUSSION 

Although, we cannot fit the 24 /tm flux measurements 
without the additional warmer component, it is clear that 
this component does not contribute significantly to the 
mass and luminosity of the region. In the regime of 
low-mass star formation, recently there have been sev- 
eral detections of very-low-luminosity objec ts (VeLLOs) 
within previously believed starless core s (lYoung et all 
12004 IDunham et all 120061 : iBourke et ail I2006D . Some 
of these objects sh ow clear outflow sign atures (L1014, 
Bou rke et all 120051 : IRAM04191+1522, IDunham et"ai] 



2006), and they are believed to be very young low-mass 



protostars or brown dwarfs likely associated with accre- 
tion disks. In the case of IR DC 18223-3, the regio n exhibits 
outflow signatures as well ([Beuther et al.l 120051 ). hence it 
is likely that the warm component in this source is pro- 
duced by an accreting low-mass protostar-disk system as 
well. An important difference between the VeLLOs and 
IRDC 18223-3 is that the luminosity of the warm compo- 
nent in IRDC 18223-3 - although it is nearly negligible 
compared to the luminosity of the cold component - is 
about 1-2 orders larger than those of the VeLLOs. Since 
we are dealing with a massive gas core that has on aver- 
age still very low temperatures and a low total luminosity, 
IRDC 18223-3 is a good candidate of being a High-Mass 
Core with an embedded low- to intermediate mass proto- 
star that is destined to become a massive star at the end 
of its formation processes. 

To constrain the status of the embedded source in 
more detail, we compare our data with the analytic and 
the radiation-hydrodynamic simulations o f the turbulent 
core model for massive star form ation (|McKee &: Tanl 
120021 120031 iKrumholz et al.ll2006l ). At early evolution- 
ary stages prior to any hydrogen or even deuterium burn- 
ing, the luminosity of the regions is completely domi- 
nated by the accretion luminosity caused by the accretion 
shocks. Th e theoretical re s ults o btained by the analytic 
approach of iMcKee fc Tanl (|2003l ) and the simulations of 
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iKrumholz et al.l (|2006f ) resemble each other well. Some 
quantitative relatively small differences are due to differ- 
ent assumptions in the initial density distributions, dif- 
ferent turbulence damping assumptions and the origin of 
the core in a shocked filament in the simulations which 
delay all result in higher accretion rates in the radiative- 
hydrodynamic simulations. Since the latter likely resemble 
real cores better, we use the simulations for our compari- 
son. 

Figure 5 in IKrumholz et al.l (|2006f ) presents the accre- 
tion rate and luminosity of the forming massive protostar 
throughout its evolution. A notable difference between the 
simulations and our observations is that the simulated lu- 
minosity is only that of the primary protostar, whereas the 
estimated luminosity of IRDC 18223-3 includes likely mul- 
tiple objects as well as potential contributions from exter- 
nal heating. However, since there is no detected UCHil re- 
gion or strong O-star within several pc from IRDC 18223-3, 
external heating contributions should be negligible. Since 
massive star formation is usually proceeding in a clustered 
mode, multiple embedded protostars below our current 
angular resolution limit (5.8" x 2.4" at 3 mm wavelength 
corre sponding to approximately 15000 AU, Beu ther et al.l 
2005) cannot be excluded. Nevertheless, it is likely that 
the luminosity of the region will be dominated by the most 
massive protos tellar object whi c h is al so confirmed by the 
simulations of IKrumholz et aTl (|2006f ). Therefore, it ap- 
pears reasonable to compare these simulations with our 
observational results. 

The luminosity obt ained for IRDC 1 8223-3 is reached 
in the simulations by IKrumholz et alj (|2006l ) already at 
extremely early times when the protostellar mass is still 
well below half a solar mass. The corresponding accretion 
rates are of the order 10 4 M Q yr -1 . Therefore, combining 
the model predictions with the observed properties of this 
region indicates that we are really dealing with a massive 
gas core that harbors an embedded, accreting low-mass 



protostar with accretion rates that are high enough that 
it will eventually form a massive star. 

Although the data do not allow us to unambiguously 
exclude that the embedded low-mass protostar may re- 
main a low-mass object potentially never forming a mas- 
sive star, the fact that the gas core is very massive, and 
that the luminosity is already relatively high (a genuine 
low-mass protostar has orders of magnitude lower accre- 
tion rates and hence a significantly lower accretion lumi- 
nosity), strongly supports the proposed scenario of a mas- 
sive star-forming region right at the onset of protostellar 
evolution and accretion. 

The example of IRDC 18223-3 shows the power of com- 
bining (sub)mm high-spatial-resolution observations with 
the recently available near-/mid/and far-infrared surveys 
of the Galactic plane by Spitzer. Combining these datasets 
allows us to constrain the SEDs of the youngest massive 
star-forming regions in unprecedented detail. This way we 
can observationally characterize the various evolutionary 
stages right at the beginning of massive star formation, in 
particular the largely unknown properties of genuine high- 
mass starless cores and high-mass cores harboring low- to 
intermediate-mass protostars like the case of IRDC 18223- 
3. We are currently performing radiative transfer model- 
ing that incorporates age estimates based on the emerging 
MIR flux but this is out of the scope of the current paper. 
For the coming years we can expect studies of larger IRDC 
samples setting constraints on the source properties in a 
statistical sense. 
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Fig. 1. — The color scales show Spit zer images at variou s wavelength. The left panel presents a three-color composite with blue 3.6 fim, 
green 4.5 (im, red 8.0 )im (adapted fromlBeuther et al,li2005T). The inlay zooms into the central core region. The middle and right panel show 
the Spitzer 24 and 70 fim images, respectively. The scaling is chosen individually to highlight the most prominent structures. Contours in 
each panel show the 93 GHz (3.2 mm) continuum emission observed with the PdBI from 1.08 mjy beam -1 (3tr) in 0.72 mjy beam -1 (2cr) steps 
HBeuther et al.ll2005l) . The axes are in R.A (J2000.0) and decl. (J2000.0). The circles in each panel present the Spitzer beam sizes and the 
ellipse in the left panel presents the PdBI 3.2 mm continuum synthesized beam. A size-ruler is also shown in the left panel. 
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Table 1 
Fluxes 



A 


S 




mJy 


24 


12.1 ±2.4 


70 


989 ± 198 


1200 


290 ± 44 


3200 


6.7 ± 1.0 




Fig. 2. — Spectral energy distribution of IRDC 18223-3. The dots with error-bars mark the detections at 24 and 70 fim on the short- 
wavelength Wien-sidc of the peak, and at 1.2 and 3.2 mm on the Rayleigh-Jeans part of the spectrum. The four stars below 10 fim show 
the Spitzer IRAC upper limits in the near-infrared. The solid line presents a two-component fit with one cold component at ~ 15 K and 
one warmer component at ~ 51 K. The two dotted lines show the two components separately. The resulting physical parameters for each 
component are labeled accordingly. 
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